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ABSTRACT 

We study the semiclassical behavior of the focusing nonlinear Schrodinger equa- 
tion in 1+1 dimensions under discontinuous "barrier" initial data and we describe 
the violent oscillations arising in terms of theta functions. The construction of 
proofs relies on 

(i) the analysis of the associated Riemann-Hilbert factorization problem 

(ii) the analysis of the resulting maximin variational problem for a Green's poten- 
tial with external field and in particular the proof of existence of a regular solution, 
which enables the construction of the so-called g-function transformation and hence 
the asymptotic deformation of the Riemann-Hilbert factorization problem to one 
that is explicitly solvable in terms of theta functions. 

In particular we show that the finite genus ansatz is generically valid for all 
times. 
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0. INTRODUCTION 

The semiclassical limit of the 1+1-dimensional, integrable nonUnear Schrodinger 
equation with cubic focusing nonlinearity has been the subject of recent investi- 
gations. Several numerical studies have appeared since 1998 ([MK], [BK], [CT], 
[C]) and rigorous analyses of the initial value problem under real analytic data 
have already appeared ( [KMM] in 2000 and [TVZ] in 2004) in cases where the 
initial data is analytic. The present paper makes use of the method and results 
in [KMM] for the study of a very particular problem with discontimious (barrier) 
data. It has been shown by A.Cohen and T.Kappeler [CK] that weak solutions 
exist for all time under barrier data and, even more, that the inverse scattering 
technique is still applicable to the integration of the problem. In fact, it is known 
([DZ2]) that there exists a unique solution in C(R+, L'^{dx) U {L°°{dx) (g) Lf^^{dt)). 
Here, we use the inverse scattering method of [CK] to pose an associated Riemann- 
Hilbert factorization problem, which we then asymptotically (as /i — > 0) reduce to 
the Riemann-Hilbert factorization problem that can be explicilty solved in terms 
of theta functions. 

The real aim of this paper is to indicate that a discontinuity in the initial data 
does not necessarily alter the behavior of the semiclassical focusing NLS problem. 
Some changes of course have to happen. For example, if the Euler system that 
appears as a formal limit of the focusing semiclassical NLS does not even admit 
a solution for small times, it is obvious that the genus zero ansatz cannot hold 
uniformly for small times. 

A natural generalization of the barrier data problem is the problem of general 
step data. For such data, the eigenvalue density p° has several nonanalyticity points 
on the "spike" where the eigenvalues accumulate to the complex plane. The analysis 
of the present paper can be immediately extended to that case. The appropriate 
contour still exists and the finite genus ansatz still applies. The proof of that fact 
is not very different from the proof of the simplest case that we present here. 
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We consider the nonlinear Schrodinger equation (l+l-dimensional, integrable, 
focusing case), on the half-hne 

(1) ^ 

u^{x,0) = uo{x), 

under barrier-hke initial data: 

uo{x) =A, -1/2 <x< 1/2, 

(2) 

= 0, otherwise. 

Here A is a fixed positive constant and ft, is a small positive constant. Eventually 
we will take h ^ 0. We will assume that h is staying away from the discrete set 
^ (2k+i)v }^ A; = 0, 1, 2, .... For simplicity, we will actually require that 



(2a) h takes values in the discrete set {- — |, fc = 0, 1, 2, ... 

KTT 



Setting 



(1) is transformed to 



H^ = -i^{u'^ui-u'u'j, 



dtp"" + d^n^ = 0, 

9*^" + 9.(^) - a,(p'')V2 = ^d.ip'^dllogpf^). 

The formal limit as /i — > is the Euler system 

dtp + d^n = 0, 

(la) /i^ 

9tM + 5x(^)-S,(p)V2 = 0. 

The initial data become p = Uq{x), ^ = 0. 

This initial value problem admits a weak solution for all time, as shown by 
A.Cohen and T.Kappeler ([CK]). Furthermore, the inverse scattering theory is still 
applicable. 
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The associated linear system is 



(3) hi;, = { J^^, ) V, 



-u*{x) iX 

where * denotes complex conjugation. Jost functions 0, ■0 are defined as column 
vector solutions of (3) satisfying the asymptotic conditions 



e" 

/ ^-i\x/h 



V' - 1 „t\x/h ) > as +00, 



Q , , OS a; — »■ — oo. 

One can define the refiection and transmission coeflacient as follows. 

Let be the transpose of {tjj2, —^Ji)- Then, following [CK] define a+, b+ by 

(I) = a+iP + b+i). 

The reflection and transmission coefficients are given by 

r+ = h+/a+, 
t+ = l/a+. 

Similarly, one can define coefficients r_ , <_ by normalizing the Jost functions at the 
opposite infinities. 
Prom [CK], we have 

(5) 

i){x,X) = ^ ^ixx/h^ , X < -1/2, 



Ae \\.L 1)1 n) (^2_,_;^2)l/2((^_i)/,j) 



, -1/2 < a; < 1/2, 



where 



(6) , {A +\)/ 

r,__ . a/. sznP^ + A^)V^(l//i)] 

^ (.l2 + ;^2)l/2 
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The coefficients a+,b+,r are given by 

a 



|a|2 + |/3|2' 



r(A) = ^- = -A- 



a cot[{A^ + A2)i/2(i//i)](A2 + A2)i/2 - iA' 
The eigenvalues are the zeros of a+. They lie on the imaginary interval [—iA, iA] 
and are given by 

A = IT], 

In other words, 

(^2- ^2^1/2 (^2 _ 2)1/2 

■5^ — = arctanl- — — + fcTr, fee Z. 

As /i — > 0, we get an asymptotic expression for the eigenvalues Afe = irjk- 
(8') {A^ - rgf/^ ~ hkiT, fc G Z. 

The limiting density of eigenvalues is then given by 

(9) 7r(^2 _ ^2)1/2- 

or, in terms of A = ir], 

(9a) p\X) = ^ 



7r(^2 + A2)l/2- 

Here the branch is chosen such that p°(A) ~ as A — > ioo. 

We note here that as a consequence of the simplifying condition (2a) we have 
/o^ p^{X)dX = —iA/n — ikh, for some integer k. Hence 

(9b) ea;p[jy^ pO(A)dA] = 1. 

This will simplify the analysis of the "local" Riemann-Hilbert problem on a cross 
centered at the origin; see section 5, paragraph 6. Another consequence of (2a) is 
that r(0) = 0. 
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It is easily seen that the associated norming constants can only take the values 
— 1, 1 (by symmetry in x) and that in fact they have to oscillate between these two 
values (by a Sturm-Liouville oscillation argument; see [KMM]). 

REMARK. The density p° blows up at iA. This is due to the local "flatness" 
and not the discontinuity of the initial data. There are some minor implications of 
this fact. The external field of the associated variational problem (see later) is still 
continuous ar iA and analytic nearby. On the other hand we do not want to allow 
that the band endpoints Xj (defined later) are equal to iA. Generically (in x, t) 
this will not happen. If, say, Aq = iA, then the anlysis of [KMM] breaks down. In 
particular it is not clear how to solve the scalar Riemann-Hilbert problem for the 
g-function. 

We end this section by defining the three Pauli matrices; we will be using them 
later. 

= { 1 o)''''={ i "o')'''' = ( -i)- 

2. THE RIEMANN-HILBERT PROBLEM 

We can now state the following Ricmann-Hilbcrt problem, following [KMM]. Let 
C be a piecewise smooth loop encircling all eigenvalues in the upper half-plane 
and lying entirely in the upper half-plane except for the point 0. Let C* be its 
conjugate, encircling the eigenvalues in the lower half-plane. Also give the following 
orientation: 

(i) the real axis is oriented from left to right, 

(ii) the loop C is oriented counterclockwise, 

(iii) the loop C* is also oriented counterclockwise. 

We use the following convention: the -|— side of an oriented contour is always to 
its left, according to the given orientation. 

We also choose C so that it approaches (from left and right) at a non-zero, 
non-straight angle with the real axis. 
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THEOREM 1 (discrete version). Let dfi = Sfe/l(^;^• — S\^) be a finite sum of 
point measures supported at the eigenvalues of the system (3), as given by (8). Let 

(10) X(A) = -(A2 + a2)1/2. 

Letting M+ and M_ denote the limits of M on S = C U C* from left and right 
respectively, we define the Riemann-Hilbert factorization problem 

M+(A) = M_(A)J(A), 

where 

2.A. + 2iA^t , ,,,,, ,AGM, 

r*(A)e H l + |r(A)|2 / 

= a2v(A*)V2, A gC*, 
limx^ooM{X) = I, 

where 
(12) 

.(A) = 



1 -i exp{j^ J log{X - r])dfi{r]))exp{-j^{2iXx + 2iXH - X{X))) \ 

1 ; 



and r is the reflection coefficient defined in (7): 

1 



(13) r(A) 



C0i[(A2 + A2)l/2(l//j)](A2 + A2)V2 _ iA ' 

Note that r(A) can have a real singularity at A = if is an odd multiple of 
7r/2. By assumption, we have excluded such values of h. 

The above Riemann-Hilbert problem admits a solution and the solution of (1) 
can be recovered from the solution of (11) as follows. 



(14) u{x, t) = 2ih limx^oo{XM'^^{X)), 

where the index 12 here denotes the (12)-entry of a matrix. 

PROOF: Standard; see Chapter 2 of [KMM]. The quantity -i exp(i(X(Afc))) 
oscillates between —1 and 1. The function —i exp{j^{X {X))) is thus an extrapolation 
of the norming constants. 
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REMARKS. 1. In [KMM] the Riemann-Hilbert problem jump involves some 
parameters denoted by K, u, J. Here we are simply choosing K = — 1,(t = 1,J = 1. 
This is compatible with the discussion in [KMM] as long as we focus our attention 
to the case a: > 0. By the obvious symmetry x — » —x this is acceptable. 

2. Obviously, the contour C can be deformed anywhere in the upper half-plane, 
as long as it passes through (we only require it to be non-tangent to either the 
real axis or the imaginary axis at 0) and does not touch the linear segment [0,2^]. 
For the asymptotic analysis, the contour C will be eventually fixed by the choice 
of the "g- function" in the next section. Similarly for its conjugate C* . In a sense 
we are seeking a steepest descent contour, on which the Riemann-Hilbert problem 
will reduce to one that is explicitly solvable. 

We note the following factorization of the jump J(A) on the real line. 




r(A)e 



-2i\x-2i\'^t 

e ^ 



= L{X)U{X), 




(15+) 



where 




A different factorization is also possible. 




r(A)e 



-2iXx — 2iX'^t 

e h 



] 



= S{X)D{X)T{X) 




where 



(15-) 





Both factorizations will be useful later. 



REMARK: From now on we will substitute = {p°{ri) + {p°)*{ri*))dri in (12), 
where p° is the asymptotic density of eigenvalues given by (9a), and supported on 
the linear segment [0, iA]. In other words, we will approximate a discrete density of 
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eigenvalues by a continuous one. This is not a trivial assumption, but it is true. We 
refer to Chapter 3 of [KMM] for a rigorous justification. Of course the justification 
of [KMM] is devised for real analytic data, but the proof goes through unaltered 
since the nice behavior of the asymptotic density of eigenvalues at the crucial point 

is not altered even when we consider steplike data. 
We then restate Theorem 1 as follows. 

THEOREM 1 (continuous version). Let = {p°{v) + {P°)* iv*))dv, where /9° 
is the asymptotic density of eigenvalues given by (9a), and supported on the linear 
segment [0,ij4]. Set 

(10) X(A) = -(A2 + A2)1/2. 

Letting M+ and M_ denote the limits of M on S = C U C* from left and right 
respectively, we define the Riemann-Hilbert factorization problem 

M+(A)=M_(A)J(A), 
where 

r*(A)e^^^n^ l + |r(A)|^ J 

= v(A),AeC, 

= (T2V(A*)V2,AgC*, 

limx^ooM{X) = I, 



where 
(120 
viX) 



1 -i exp{\ J log{\ - ri)dn{ri))exp{-\{2i\x + 2i\H - X{\))) \ 

1 ; 



(13) r(A) = -~A- 



COt[{A^ + A2)l/2(l//i)](A2 + p)l/2 _ 

The above Riemann-Hilbert problem admits a solution and the solution of (1) 
can be asymptotically recovered from the solution of (11) as follows. As — > 0, 

(14') u(x, t) ~ 2ih ZimA^<x)(AM^2(A)), 
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where the index 12 here denotes the (12)-entry of a matrix. 

PROOF: See Chapter 3 of [KMM] in the case where the maximizing contour of 
the variational problem (see Appendix 2) does not touch the spike [0,iA]. For the 
general case, follow the discussion in Appendix 2. 

3. ASYMPTOTIC ANALYSIS OF THE RIEMANN-HILBERT PROBLEM. 

The ideas underlying the discussion of this section are the following. 

(i) Our given Ricmann-Hilbcrt problem needs to be asymptotiealy deformed to 
an explicitly solvable one. To do this, some analyticity properties of the jump matrix 
are crucial. Proofs of the passage from one Riemann-Hilbert problem to another 
rely heavily on the equivalence of Riemann-Hilbert problems and special singular 
integral equations (see [DZl] for the first rigorous description of the deformation 
method). 

(ii) One of the crucial deformations involves a certain conjugation and the intro- 
duction of an appropriate "g- function" (first introduced in [DVZl]; see also [DVZ2]). 
See (17) below. 

(iii) For the focusing nonlinear Schrodinger problem and because the associated 
Lax operator is not self — adjoint the appropriate introduction of the g-function 
requires finding a "steepest descent contour"; see (18) below. The term "nonlinear 
steepest descent method" often applied to the Riemann-Hilbert deformation theory 
initiated in [DZl] thus acquires full meaning. For this we follow the ideas of [KMM], 
[KR]; see also Appendix A. 2. 

We begin with the following observation. Consider the reflection coefiicient given 
by (13). For A in the upper half-plane, at least away from the real line and the 
eigenvalues given by (8), we have 



while for A in the lower half-plane, at least away from the real line and the eigen- 



(16) 



r(A) 



-iA 



(A2 + A2)l/2+A' 
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values given by (8), we have 
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(16*) 



r(A) 



iA 



(^2 + A2)l/2_A" 



This means that r(A)e 




is exponentially decaying (growing) in the upper 



half-plane, at least away from the real and imaginary axes, as long as Re(X) < —x/2t 



the lower half-plane, at least away from the real and imaginary axes, as long as 
Re{X) < -x/2t {Re{X) > -x/2t). 

The above suggests that using the factorizations of the Riemann-Hilbert problem 
defined by (15-), (15-I-) and applying the obviously suggested deformations, the 
jump across the real line should be reduced to a diagonal matrix independent of 

X, t. 

To do this of course, we must ensure that the solution of the Riemann-Hilbert 
problem is uniformly bounded (in x, t) as h ^ 0. This is actually not true as will be 
seen later. What is true however, is that a judicious conjugation of the jump matrix 
will "deform" it to a new Riemann-Hilbert problem whose solution is uniformly 
bounded (in x, i) as /i — >■ 0. At that point we can neglect all terms involving r. 
Furthermore, the very same conjugation will ensure that the new Riemann-Hilbert 
is explicitly solvable (asymptotically). The deformed Riemann-Hilbert problem will 
be defined with respect to a steepest descent contour C. 

More precisely, we introduce the change of variables 



(i?e(A) > -.T/2t), while r*(A*)e 



2iXx + 2i\'^t 



is exponentially decaying (growing) in 



h 



(17) 



Q(A) =M(A)e^, 



where (73 = diag{l, —1) (a Pauli matrix) and the complex-valued function g is 
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constrained by the following conditions: 

g{X) is independent of h. 
g{X) is analytic forX G C \ (C U C*). 

g{X) — > as A ^ 00. 

(17a) 

g{X) assumes continuous boundary values from both sides of CuC*, 

denoted by g+{g-) on the left (right) of C U C* . 

g{X*) + g{Xy = for all XeC\{CU C*). 

These conditions of course do not define g uniquely. They will be augmented by two 
conditons below, that will also implicitly define an admissible steepest descent contour C. 

The assumptions above permit us to write g in terms of a measure p defined on 
the contour C. Indeed 

(18) g{X) = log{X - v)piv)dv, 

^ J CISC" 

for an appropriate definiton of the logarithm branch (and for a; > 0; if x < there 
is a sign change in [KMM] but of course we can restrict ourselves to the case a; > 
here because of the symmetry of NLS). 
For A G C, define the functions 

^(A) ■.= i{g+{X)-g.{X)\ 
</.(A) := f log{X - rj)p\r]) drj+ [ log{X - 7?)p°(r?*)* drj 

Jo J-iA 



(18a) 

+2iXx + 2iX^t + in / p"{r]) dr] - g+{X) - g-{X). 



A 

Now we spell out the two important conditions which determine the (steepest 
descent) contour C (not uniquely) and the function g. 

p{ri)dr} is a real measure, 

(18b) 

5R(</.(A)) < 0, 

i.e. a "measure reality" condition and what can be interpreted (see Chapter 8 of 
[KMM]) as a "variational inequality" condition. In fact, one can eventually show 
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that the measure p{rj)dr} has to be nonpositive: strictly zero in the "gaps" and 
strictly negative in the "bands"; see section 4 later. 



REMARKS. 1. The fact that there exist a contour C and a function g defined 
implicitly by (17), (17a), (18), (18a), (18b) is highly non-trivial. In [KMM] we had 
not proved this fact for all times. Instead it was only proved for small (finite) times 
and it was simply assumed for larger times. A result for larger times appeared only 
later in [KR] (see Appendix 2 for a discussion). Having the results of [KR] (or more 
accurately a modification of them as in Appendix 2) we do not need to dwell too 
much in the details of the construction of C, g any more. 

2. The function g and the steepest descent contour C depend only on the density 
of the eigenvalues, They are independent of the reflection coefficient r. 

By (17) the Riemann-Hilbert problem for Q is Q+(A) = Q-{\)vq{\), where 



i exp{^ J log{X - r])diJ.{r]))exp{^{-g+ - g- + 1i\x + ^i^^t + X*(A*))) 

/ 



Also lim\^ooQ{X) = I- 

Now, we will be able to treat the terms involving the reflection coefficient r 
by arguing as follows. First we note that g is purely imaginary for real A; this 
follows from (18)-(18b). So, write g = iip, so that V(A) e M for A e M. Let 
C(a;, t, X)='tjj- \x - X^t. Clearly C(A) G M for A e R. 

Next, divide the real line into a union of (finitely many) intervals, say Jfe, siich 
that for A € interior(Jk), either ^ > or ^ < 0. Denote by J,^ the intervals in 
which ^ > and by Jf the intervals in which ^ < 0. Naturally, the (nonzero) 
endpoints of the intervals Jk are given by the condition ^ = 0- Note that C is 
real analytic in M \ 0, and C ~ —X'^t, as A — > ±oo. The finiteness of the number of 





g- -9+ 

e h 




e 



g+-g- 

h 
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intervals follows from the analyticity of (, at least away from 0. But also near 0, 
we will show that ^ <0. 

We first show that when A = 0, then = — oo. Recall that 



aW = \ I log{X - r])p{r])dr], 



/cue* 

where CUC* is oriented counterclockwise and p(r\)dr\ is a non-positive measure on 
C which is strictly negative in a non-trivial subset of C and is extended to C* by 
the condition 

P^rf) = {p{v)r. 



Differentiating, one gets 

dgjX) 
dX 

and at A = 0, 



^ Jcuc 



dg{X) -1 f p{t]) 



If P(JI 
' Jcuc V 



dr}. 



dX 2 

Making use of the symmetry with respect to complex reflection and remembering 
that the orientation is counterclockwise for both C and C* one gets 



The integrand is strictly negative and the integral diverges because p is nonzero at 
zero. It follows that at A = 0, = —ioo, so ^ = — oo. 
Now, if A is real and close but not equal to 0, 

dX Jc \X - T]]^ 

so ^ < 0. Hence the point A = belongs to the interior of some interval Jj~. 

Consider first the intervals J^^. By the Cauchy-Ricmann relations we have 
d{Im()/d{ImX) > 0, for A G Jj^ , in the positive direction perpendicular to the 
real line. This means that in an area of the upper half-plane, close to Jj^ , the real 
part of i( is negative, so exp{i(/h) is exponentially decaying. Similarly, in an area 
of the upper half-plane, close to Jf , the real part of i( is positive, so exp{—i(/h) 
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is exponentially decaying. In the particular case of A = it is also easy to check 
that exp{—iC,/h) is exponentially decaying as A moves upwards along the positive 
imaginary axis. 

Let us now introduce the following lens-like contours. For each interval J;"^ con- 
sider small piecewise linear deformations upwards, say J^"^'"^, keeping the end points 
fixed, but otherwise lying entirely in the upper half-plane. Similarly, consider small 
piecewise linear deformations downwards, say J'^''^^ , keeping the end points fixed, 
but otherwise lying entirely in the lower half-plane. For each interval J^^ consider 
small piecewise linear deformations upwards, say J; keeping the end points 
fixed, but otherwise lying entirely in the upper half-plane. Similarly, consider small 
piecewise linear deformations downwards, say J^"'*^", keeping the end points fixed, 
but otherwise lying entirely in the lower half-plane. All orientations are compatible 
with Re{\) increasing. We also make sure that J^'"^, J;"^'*^", Jp^^ , J^"'"^" cut the 
real line at angles ^ 0, 7r/2. See Figure 1. 




IR 



Fig.l. Part of the extended contour, with lenses 

Note here that one of the lens-like contours J;~'"^ will cut the contour C while 
one of the lens- like contours J;"'''" will cut the contour C* . 

Let the open regions bounded by and J^'^^ be denoted by D^'^^ and the 
open regions bounded by and J^''^"' be denoted by Df'''''^. Similarly, let the 
open regions bounded by J;~ and J;~'"^ be denoted by D^'^^ and the open regions 
bounded by Jf and J^"'**" be denoted by D;"'"^". We make sure that Re{iQ < 0, 
in U J;+'"^ U £>("''^" U J;"'"^". Similarly, we make sure that Re{iC) > 0, in 
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£>;-'"^ U J;-'"^ U U 7+'"^". 

Next, we can make use of the factorizations given in (15+), (15-). Set 

Z{X) = g(A), for XeC\ U/[5+'"^' U 5+'''" U L)-^"" U 

Z{X) = Q{X)U-\X), for X € UiD+'^'p, 
(19) Z{X) = Q{X)L{X), for X e U£>+''^", 

Z{X) = Q(A)r-i(A), for X G U;i^r'"^ 
Z(A) = Q(A)5(A), /or A e UZ);"''^". 

The Riemann-Hilbert problem for Z is described by 
(20) 

Z+(A) = Z_(A)i;z(A), 
J -^(^)^ ^ " j , /or A G U J+'"^ 

z;z(A) = ^ r-*(A*)e"'^"^"fe"^''^'' l) ' •^''^ ^ ^7;+'"^", 

I'Z(A) = ( r''(A'') -2g + 2iXx + 2i\^t j , /or A e U; J;~ 

V l+|r(A)pe 

vz{X) = i I TWF^ ^ " ], forXe UiJr"'', 
e^^-TT^ -iexp{\jlog{X-ri)d^i{vi))exp{-\{-g+-g-+2iXx + 2iXH-X{X)))\ 



vz{\) 



vz{X) = 



9- -9+ 

e—fr^ 



for A e C, 



9-L -9- 



i exp{=^ J log{X - 'n)dii{ri))exp{\{-g+ - g- + 2iXx + 2iXH + X*{X*))) e— 

for XgC 

<'-ir''"' ,,,%,). A. u,, 

with the normaUzation Um\^aoZ{X) = I. For a part of the extended contour, away 
from 0, and between the points ai and a(+i, including the small lenses, see Figure 
1. 

The idea for the next transformation is to get rid of the jump across the real 
line. 

We define d as follows. Let d be analytic in C \ U J^" , such that the normal limits 
from above and below the real line d± exist (apart from the endpoints) and such 



I ' 
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that 

J+(A) = d-(A)(l + |r(A)p), A e Uiinterior{J,^), 

lim\^ood{\) = 1. 

One can even give an explicit formula for d. Denote = [ai, bi]. Then 

(21) d{X) = exp [121 log{l + \r{C)\')^-^-^], A e C \ UJf 

where the principal branch of the logarithm is chosen. 

It is clear that the definitions above are not singular and indeed d,d~^ are 
bounded in the complex Riemann sphere. In fact, it is easy to check that both 
d{X) and [(i(A*)*]~^ satisfy the conditions of the scalar Riemann-Hilbert problem. 
Hence they must be equal, from which it easily follows that 

|d_p(l + |r(A)n = l, forXeUiJr, 

(21a) \d-\ < 1, for A € UjJf, 

< 1 + sups,\r{X)f < 00, for X G U; Jf, 

and hence, by the maximum principle, |d(A)| < 1 + sup]^\r{X)\'^ < oo. So \d{X)\ is 
bounded uniformly in the complex plane. Similarly, |d~^(A)| is uniformly bounded 
in the complex plane. 

Near the points a; the local behavior of d off the jump contour is 

(21b) d r-. di{X - aif' , 

and near bi, 

(21c) d^ fiiX-biY^', 

where ui = ^(1 + log\r{ai)\'^), jj^i = ^(1 + log\r{bi)\'^) and di,fi are independent 

of A. In fact 

(21d) 

. ,£ - ..^-(M , ,,,, £■ , Ko,^,^,, 
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as can be shown by integration by parts. The next transformation is then 

Y{X) = Z{X)A{X), 

where 



A(A) 



d{X) 
J- 1(A) 



(22) 

Then Y_^_ = Y_vy, where 
(23) 

/ e'^^T^ -id-'^ exp[Ujlog{X-r])dn{r])+g-+g--2iXx-2iXH + X{X))]\ 

VY{X) = g_-g , 

V e—nr^ J 

for A e C, 



vy{X) = 



9+ —a— 

e-^T— 0_ 
i exp[\{- / log{X - ri)dii{vi) - g+ - g_ + 2iXx + 2iXH + X*{X*))] 



for X e C*, 



1 -r{X)d e^ j , /or A € U, J+'"^ 
vy{X) = ^ ^^^^^^~_2 j^2g+2ixx+2ix^t , for Xg U/J^''^", 



'"y{X) = J-2 r*(A*) -23+2.Ax+2.A''t , /or A G U;J, 

^ ^ l+|r(A)|2'^ ' 

2g-2iXx-2i\'^t 



up 

I ' 



vy{X) = f J T^RAF^^^ ' j , /or A G Jf''^". 



Also /imA^ooi^(A) = /. 

At this point, we still have a Riemann-Hilbert that is equivalent to the original 
one (11)-(12), at least accepting the discrete-to-continuous passage in Theorem 1. 
Now, we can finally start considering the limit /i — > 0. Indeed, all terms involving 
r(A) and not supported on the real line can be neglected, not because r itself is 
small (it is not, see (16)-(16*)) but only because it always appears multiplied by 
something exponentially small. 

In general jump matrices of the form / + exponentially small can be neglected 
asymptotically as long as it is proved that the solution is uniformly (in x, t) bounded 
as /i — + 0. We shall see eventually that this is the case. 

Assuming for the moment that this is true, we can simply delete the non-real 
part of the contour, at least away from the endpoints of the intervals Jk- Eventually 
(see section 5 later) we can also delete the remaining small crosses centered at such 
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points. We end up with a matrix valued function W, such that W r^Y near infinity, 
and 



(24) W+{X)=W-{X)vw{X), 



where 
(25) 

vwW = 

Vw{X) 



e ^ —id exp{ i / log{X - ri)diJi{r]))exp{-\{-g+ - §- + 2iXx + 2iXH - X{X))) \ 
e—^ J 

for A G C, 



3J- -9- 

e^^T— 



i exp{^ J log{X - 'n)dn{rf))exp{\{-g+ - g- + 2iXx + 2iXH + X*{X*))) e—^ 

for X € C*. 

Also limx^ooW{X) = I. 

We will eventually check that the solution of the Riemann-Hilbert problem for 
W exists and is uniformly (in x, t) bounded as /i — > 0. This will justiiy neglecting 
the exponentially small terms in (23). The passage from (23) to (25) will then be 
a posteriori justified, again at least except for the remaining small crosses centered 
at a finite number of real points. Concerning the crosses see section 5 (Remark 5) 
and Appendix A.l. 

After applying the above transformation Y — > W, we must use another lens 
transformation to simplify the jump across the conotur CUC*. We will not describe 

this new lens transformation, since the discussion is exactly as in Chapter 4 of 
[KMM]. We will simply state the end result of these lens transformations and the 
transformation Y W. The result is the so-called outer problem. We shall show 
in the next section how it can be treated along the lines of Chapter 4 of [KMM]. 

THE OUTER PROBLEM. 

We first define the analytic arcs Ij,I*, j = l,...,G/2 as follows (they come in 
conjugate pairs). Let the points Xj, j = 0, ...,G, in the open uper half-plane be 
the branch points of the function g. (The fact that there are G + 1 of them is 
a consequence of the definition of g, according to the "finite genus ansatz"; see 
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below.) All such points lie on the contour C and we order them as Aq, Ai, Ag, 
according to the direction given to C. The points Aq, A^, Aq are their complex 
conjugates. Then let /q = [0, Aq] be the subarc of C joining points and Aq. 
Similarly, Ij = [\2j-i,^2j], j = l,..-,G/2. The points Xj, j = 0, ...,G lie in the 
open upcr half-plane and they are determined by a set of transcendental equations 
that follow directly from the definitons of g and C (cf. Remark after Lemma 5.1.5 
of [KMM]). The connected components of the set C \ U /*) arc the so-called 

"gaps" , for example the gap Fi joins Aq to Ai, etc. The subarcs Ij are the "bands" . 

ASSUMPTION (G). For simplicity, we will make the generic (in x, t) assumption 
that Xj i^iA, j = 0, ,G. 

The finite genus ansatz implies that for each x, t there is a finite positive integer 
G such that the contour G can be divided into bands and gaps as above. In 
fact, it follows from the conditions defining p, G that the measure reality condition 
{p{r})dr] real; see (18b)) splits into a measure strict negativity condition in the 
bands and a measure zero condition in the gaps. Furthermore, the function 9{X) of 
(18a) defined on C is constant on each of the gaps Fj, taking a value which we will 
denote by 6j, while the function (j) of (18a) is constant on each of the bands, taking 
the value denoted by aj on the band Ij. For the justification of the finite ansatz 
under the barrier data see Appendix 2. 

We are seeking a matrix O, which is analytic everywhere except across the con- 
tour C\Fcr/2-i-i and its conjugate, with limits that are i'^(C\F(j/2-M)) converging 
to the identity at infinity and such that 

n (\\-n (\\( ^ id-'^exp{-iak/h)\ ^ ^ r, . r* 

0+(A) - 0_(A) ^ ^J2g^p(^^^/^) J , A e 4 U 4, 

(26) , ^ = 0,1,...,.A 

0.i^) = 0.iX,{ -^f/^^ e..(-V))'^^^^^^^' 

k^l,...,G/2. 

We here recapitulate the sequence of matrix deformations inroduced so far: 



M{discrete) -»• M{continuous) -^Q^Z^Y^W^O. 
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The first problem in the sequence, for M (in its discrete version), is equivalent to 
the inverse scattering problem for NLS. The last Riemann-Hilbert problem, for O, 

(26) , will be solved explicitly via theta functions (see (33) of section 5). 

4. THE RESULT 

As explained in Appendix 2, the finite genus ansatz holds for the semiclassical 
asymptotics under barrier initial data, at least generically. 

Assuming the finite genus ansatz, the x, t-plane can be divided into (possibly 

empty) open regions Rg, G = 0,2, 4, (here G/2 + 1 is the number of components 

of the support of the equilibrium measure of C, see Appendix A. 2), together with 
their boundaries, such that within each region the asymptotics of the solution of 
(1) with barrier data as in (2) can be given as follows. 

THEOREM 2. Let xq, to lie in region Rq. The solution u'^{x, t) is asymptotically 
described (locally) as a slowly modulated G+1 phase wavetrain. Setting x = xo+hx 
and t = to + hi, so that xo,to are "slow" variables while x,t are "fast" variables, 
there exist parameters 

a,U = {Uo,Ui,....,UgY' , k = {ko,ki, ,^0)^, w = {wo,wi, ,wg)'^, Y = 

{Yo,Yi, ,Yg)'^, Z = {Zo,Z\, ,ZgY^ depending on the slow variables xo 

and to and possibly h (but not x, i) such that 

u^{x,t) = u^'ixo + hx,to + hi) - a(a;o,to)e'^°('''"*°)/''e*(*^°(^'"*°)*-'"°(^°'*°)*' 

(27) ^2^^ ~^ Q{Y{xo,to) + iU{xo,to)/h + i{k{xo,to)x - w{xo,to)i)) 

Q{Z{xo, to) + iU{xo,to)/h + i{k{xo, to)x - w{xo,to)i)) 
Unlike the analogous formula in [KMM], we allow here a dependence of the param- 
eters on h. But of course, we can always rearrange terms to arrive at a formula like 
(27) where the parameters are independent of h. 

All parameters can be defined in terms of an mider lying Ricmann surface X. 
The moduli of X are given by Xj, j — 0, ....,G and their complex conjugates 
A^, 2 = 0, G. The genus of X is G. The moduli of X vary slowly with x,t, i.e. 
they depend on xo,to but not x,t. For the exact formulae for the parameters as 
well as the definition of the theta functions we present the following construction. 
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The Riemann surface X is constructed by cutting two copies of the complex 
sphere along the slits loD lQ,Ij,Ij,j = l....,G, and pasting the "top" copy to the 
"bottom" copy along these very slits. 

In exact analogy with the discussion of [KMM] (as exemplified in Figure 4.5 
on p. 62 of [KMM]) wc define the homology cycles aj^bj, j = 1,...,G as follows. 
Cycle ai goes around the slit Jg U /g joining Aq to Aq, remaining on the top sheet, 
oriented counterclockwise, 02 goes through the slits I-i and Ii starting from the 
top sheet, also oriented counterclockwise, as goes around the slits /-i,/o U 
remaining on the top sheet, oriented counterclockwise, etc. Cycle bi goes through 
lo and Ii oriented counterclockwise, cycle 62 goes through 7_i and /i, also oriented 
counterclockwise, cycle 63 goes through 7_i and I2, and around the slits /_i,/o U 
/q , /i , oriented counterclockwise, etc. 

Let 

G 

R{X)' = l[{X-Xk){X-Xl), 

fe=0 

choosing the particular branch that is cut along the bands Ii~ and and such that 

On X there is a complex G-dimensional linear space of holomorphic differentials, 
with basis elements t'fe(-P) for k = 1, . . . , G that can be written in the form 

G-i 



where Rx{P) is a "lifting" of the function i?(A) from the cut plane to X: if P is 
on the first sheet of X then Rx{P) = R{X{P)) and if P is on the second sheet of 
X then Rx{P) = — i?(A(P)). The coefficients Ckj are uniquely determined by the 
constraint that the differentials satisfy the normalization conditions: 



(b Vk{P) = 271-16 jk- 
J ai 
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Prom the normalized differentials, one defines aGxG matrix H (the period matrix) 
by the formula: 




It is a consequence of the standard theory of Riemann surfaces that H is a sym- 
metric matrix whose real part is negative definite. 
In particular, we can define the theta function 

<d{w) := exp{—n^Hn + n^w), 

where H is the period matrix associated to X. Since the real part of H is negative 
definite, the scries converges. 

Wc arbitrarily fix a base point Pq on X. The Abel map A : X ^ Jac{X) is then 
defined componentwise as follows: 

Ak{P;Po) := / MP'), k=l,...,G, 

JPo 

where P' is an integration variable. 

A particularly important element of the Jacobian is the Riemann constant vector 
K which is defined, modulo the lattice A, componentwise by 

where the index k varies between 1 and G. 

Next, we will need to define a certain meromorphic differential on X. Let fl{P) 
be holomorphic away from the points ooi and 002, where it has the behavior 

0(P) = dp{X{P)) + , P - 001, 

n{p) = -dp{x{p)) + o (^) , P - 002, 

and made unique by the normalization conditions 

/ n{P) = o,j = i,...,G. 
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Here p is some polynomial. In the present context see section 5, equation (31) for 
its definition. 

Let the vector U e C*^ be defined componentwise by 

Jb, 

Note that has no residues. 

Let the vectors V\,V2 be defined componentwise by 

Vi,k = {Ak{K+) + + Ak{Xl+) + ■■■ + Ak{XG+)) + ^fc(oo) + 7,1 + 

V2,k = {Ak{\\+) + Ak{\2+) + Ak{Xl+) + ■■■ + Ak{XG+)) - Ak{<x>) + m + 

where k = I, ...,G, and the + index means that the integral for A is to be taken on 
the first sheet of X, with base point Aq. 
Finally, let 

kn = dxU„, Wn = -dtUn, n = 0, . . . , G, 

where 

Y = -A{oo) -Vi, Z = A{oo) - Vi, 

and Uo = —{61+ chq) where 61 is the (constant in A) value of the function 6 in the 
gap Fi and ao is the (constant) value of the function cp in the band Iq- The fact 
that these values are actually constants in A follows from the conditions defining g 
and C. 

Now, the parameters appearing in formula (27) arc completely described. 
We simply note here that the Ui and hence the ki and Wi are real modulo 0{h). 
We also note that the denominator in (27) never vanishes (for any xo,tQ,x,i). 

REMARK: Because C depends only on the eigenvalue density and not on 
the reflection coefficient r, the constructions of C, the Riemann surface X, the 
holomorphic differential VkiP), the Abel map A and the theta functions are all 
independent of r. The only contribution of r comes through the factors (P,d~'^ 
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in (26) and the factor (r{Xo) in (27). This is why our discussion in section 4 is 
virtually repeating verbatim the analogous discussion of section 4 in [KMM] . 

5. REMARKS AND PROOF OF THEOREM 2. 

1. Formula (27) is locally a so-called finite gap expression. It describes violent 
oscillations of bounded amplitude but high frequency 0{l/h). 

2. As in [KMM] weak limits of densities exist: p = limfi^o\u^\'^ and /i = 
limh^o—Y^{u^u^ — u^u^). Those limits are actually strong in the genus zero region. 

3. Naturally, since the initial data is discontinuous and the limiting Euler system 
(la) is elliptic (see [KMM]) one expects the break-time of the limiting system to be 
zero. This is indeed the CcLSGj cLS numerical experiments by H.Ceniceros and F.Tian 
have shown [T], or as can be shown analytically by considering the limiting Euler 
system directly: from the second equation of (la) it is obvious that fit is infinite 
ai X = ±^,t = 0. Of course a genus zero region still exists but the first caustic 
(the boundary between the genus zero region and higher genus regions) touches the 
t = axis of the x, t-plane a.t x = ±1/2. 

4. The proof of the results in [KMM] makes use of the assumption that the 
eigenvalue density can be analytically extended in the upper half-plane with the 
spike where eigenvalues accumulate deleted. This is the case here, see formula 
(9a). The branch root singularity at iA is integrable. Integrals like J p^dr] or 
/ log{X — r])p^{ri)dr} can still be deformed and the Cauchy theorem holds. On the 
other hand, we assume that x, t are such that ^ iA, j = 0, 1, ....g, a generic case. 

5. A priori, neglecting the exponentially terms in (23) only allows us to "delete" 
arcs J;^'"^, J;^'''". J; ' ; J; ' ouly away from thc Tcal cndpolnts of J;^ . J;^ . For en- 
tirely rigorous justifications of the deletion of the small crosses centered at each such 
point remaining after the deletion of the bulk of the arcs J;"*"'"^, J^''^"', Ji~'^^, J^''^^ , 
one must construct local parametrices of the Riemann-Hilbert problem and make 
sure they match with the solution of (25) away from the endpoints. This is a 
procedure that is by now standard in the literature; more details are given in the 
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Appendix. In fact, essentially the same situation has appeared in [DZl]. The "lo- 
cal" Riemann-Hilbert problem can be solved via parabolic cylinder functions. (In 
[DZl], of course, the contribution was not negligible, since one was trying to eval- 
uate the term of order 0{t~^^^) of the asymptotics. Also in [DZl] the reflection 

coefficient is independent of the parameter 1/t going to zero, but this is irrelevant 
since the local Riemann-Hilbert problem is solved exactly via parabolic cylinder 
functions.) 

6. As in [KMM], one needs to provide a local parametrix near the origin and then 
match it with the solution of the "outer problem". This can still be done, using 
the Predholm theory described in [KMM]. The "cyclic" relation, that the product 
of limits of jumps at is the identity, still holds. In view of (9b) the cyclic relation 
follows easily from the analogous relation in the reflectionless case. The discussion 
in sections (4.4.3), (4.5.1) and (4.5.2) of [KMM] can be then followed verbatim. 

7. The Bohr-Sommerfeld condition (8) is not quite the same as the condition 
postulated in [KMM] defining the so-called soliton ensemble, which translates as 

for our present problem (l)-(2). Now, one can notice that the difference between the 
two conditions gives rise to a uniform error of higher order 0{h^). The analysis in 
[KMM] (Chapter 3) of estimates needed for the passage from a "discrete" Riemann- 
Hilbert problem to a " continuum" Riemann-Hilbert problem is not altered by this 
innocent modification. 

8. Prom formulae (6)-(7) it is obvious that a(0) = if Ah is an odd multiple of 
7r/2. So, for a particular sequence of h going to zero we have a spectral singularity 
at the real point 0. The Riemann-Hilbert problem jump becomes singular at 0. By 
assumption we have excluded such values of h. We plan to study the effect of a 
real spectral singularity on the semiclassical behavior of the focusing NLS in a later 
publication. 

9. The proof of Theorem 2 now follows the discussion of Chapter 4 of [KMM]. 
One minor change is the extra factor appearing in the off-diagonal terms of (26). 
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This factor can be taken care of by the auxihary scalar Riemann-Hilbert problem 
(4.38) in section 4.3.1 of [KMM]. The proof goes through with only a minor change: 
ak has to be substituted by ak + 2ih logd. Of course d has to appear in the final 
formula (27). Granted that the term hd is no more constant on bands, but then 
the auxiliary scalar Riemann-Hilbert problem is still explicitly solvable. 
More precisely, we introduce the scalar problem: 

H+{X) - H_{X) = -Ok, A e Tfe U fc = 1, . . . , G/2 , 

(28) 

H+{X) + i?_(A) = -{ak + 2ihlogd), A e 7fc U 7^, fc = 0, . . . , G/2 , 
not specifying any condition at infinity yet. Consider the matrix defined by 

P(A) = 0(A) exp{iH{X)a3/h), 

where O is the solution of the outer Riemann-Hilbert problem (26). It is straight- 
forward to verify that the matrix P{X) has the identity matrix as the jump matrix 

in all gaps Tk and F^. Since the boundary values of 0(A) and H{X) are continuous, 
it follows that P{X) is in fact analytic in the gaps. In the bands, the jump relation 
becomes simply 

P+(A)=iP_(A)ai, 

so the jump relation is the same in all bands. Next, suppose that (3{X) is a scalar 
function analytic in the A-plane except at the bands, where it satisfies P+{X) = 
—i[3-{X). Suppose further for the sake of concreteness that /3(A) ^ 1 as A ^ oo. 
Then, setting 

(29) V{X) = /3(A)P(A), 

we see that the jump relations for V{X) take on the elementary form: 

(30) V+{X) = V- (A)ai, A G Ufc(7fe U 7^*). 



Our purpose in reducing the jump relations to this universal constant form is that 
it can be explicitly solved in terms of theta functions. 
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But let US describe the scalar functions H{X) and /3(A). We get 



G/2 



A — An -n- A — A2fe-i A — Aj^ 



A - Ao ^^-fj A - A2fe_i A - A2fe 

and for /3(A) we select the branch that tends to unity for large A and that is cut 
along the bands Ik and I^. It is easily checked that /3(A) as defined here is the only 
function satisfying the required jump condition and normahzation at infinity that 
has continuous boundary values (except at half of the endpoints). To find H{X), 
we introduce the function i?(A) defined by 

R{Xf = f[{X-Xk){X-Xt), 

fe=0 

choosing the particular branch that is cut along the bands Ik and I^ and satisfies 

hm = -1, 

This defines a real function, i.e. one that satisfies i?(A*) = R{X)* . At the bands, 
we have R+{X) = — i?_(A), while R{X) is analytic in the gaps. Set 

H{X) = k{X)R{X), 

where 

^^"^ dr] 1 f (ttn + 2ihlogd{r])) drj 



^^^^ - L 5/"/r„ur, (A -v)R{ri) ^ L ^Ji^.i, 
We see that A;(A) satisfies the jump relations: 



(A-r,)i?+(r7) 



fc+(A)-fc_(A) = --^, Aer„ur; 

u (\\ u (\\ Qn + 2ihlogd{X) 

k+{X) - fc_(A) = R^) ' A e in U i„, 

and is otherwise analytic. So H satisfies (28). 

The function k blows up like (A — A„)~^/^ near each endpoint, has continuous 
boundary values in between the endpoints, and vanishes like 1/A for large A. It is 
the only such solution of the jump relations. The factor of R{X) renormalizes the 
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singularities at the endpoints, so that, as desired, the boundary values of H{X) are 
bounded continuous functions. Near infinity, there is the asymptotic expansion: 
H{X) = HgX^ + Hg-iX^-^ + --- + HiX + Ho + O(A-i) 

(31) 

= p(A) + 0(A-i), 

where all coefficients Hj of the polynomial p{X) can be found explicitly by expanding 
i?(A) and the Cauchy integral fc(A) for large A. It is easy to see from the reality of 
6j and aj that p{X) is a polynomial with coefficients which are real modulo 0{h). 
So the matrix function V defined in (29) has the following asymptotics at infinity: 

(32) V{X)exp[-ip{X)a3/h] = 1 + 0{X-^). 

This together with the jump relations (30) defines a Riemann-Hilbert problem for 
V that can be explicitly solved via theta functions. Equivalently P and O can be 
explicitly expressed in terms of theta functions. For example, it is now elementary 
to check that the solution of the outer Riemann-Hilbert problem (26) for O is given 
by the formulae: 

6- (A) e(A(oo) - Vi) e{A{x) -Vi + iu/n) 



On (A) 



/J(A) e(A(A) - Vi) Q{A{oo) - Vi + iU/h) ' 



n (\) = ^^(^) ,2iR(x)kM/n e(^(^) - ^i) e(-A(A) -v, + tu/h) 
''^ ' m &{-A{x) - VI) e(^(^) -v, + iu/h) ' 

&+(A) „_2.i^(A)fe+(A„)A0(-^(oo) - 1^2) e(/A(A) - ^2 + iC//?l) 



021(A) = ^T^e 

022(A) = 



/?(A) Q{A{X)~V2) Q{-A{oo)-V2+iU/hy 

b-{X) e{-A{oo) - V2) Q{-A{X) -V2+ iU/h) 



where 

(34) 6±(A) 



/3(A) e(-^(A) - V2) e{-A{oo) -V2 + lU/h) ' 

i?(A)±(A-Ag)(A-Ai)...(A-A^) 



2i?(A) 

Prom the explicit solution, using formula (14) with the obvious substitution of O 
for M one derives formula (27). This completes the proof of Theorem 2. 

We end this section by once more recapitulating the sequence of matrix defor- 
mations used: 



(35) M{discrete) M{continuous) ^Q^Z^Y^W^O. 
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The first problem in the sequence, for M (in its discrete version), is equivalent to 
the inverse scattering problem for NLS. The last Riemann-Hilbert problem, for O, 
(26), was solved explicitly via theta functions. Since, as is seen, the solution O of 
the outer problem is uniformly bounded in x,t, with limits 0+, 0_, as h ^ 0, 

so arc W and Y and the (asymptotically valid) transformations from y to O and 
back are a posteriori justified. 
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APPENDIX 1. THE CROSS PROBLEM. 

The transformation from the Riemann-Hilbert problem (23) to the problem (24)- 
(25) requires two steps. First, the deletion of the lens contours away from the points 
ai,bi. This is immediate because the jump matrices are uniformly exponentially 
small perturbations of the identity. 

Second, one needs to consider the small remaining crosses centered at the points 
ai, hi (see the remark of section 5). Since the jump matrices are not uniformly small 
there, one needs to find "local" parametrices. In other words, one needs to solve 
the local Riemann-Hilbert problems. 

For example, after translation all problems centered at ai look as follows: to find 
a matrix L which is analytic in C \ F where F is the cross shown in Figure A.l, 
centered at 0. The actual angles between the four half-lines emanating from are 
not important as long as every half-line is in a different quadrant. 




Fig.A.l. The Riemann-Hilbert problem on a cross centered at zero. 
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The jumps for L are 



Z/+ = L_?;l(A), where 



(A.1) ^4A)=( ^-2_^;,-24 -"'"^^-"^'^^' J),/orAeJ2, 

I -dh>^"^"e^^ ),/orAeJ3, 

-"lW = ^-2^*^-2L.,ig^M±Hi^5±2i^ , /or A e J4. 

Also lim.\^ooL{X) = I. 
Here 

91 = ^(A = a;)>n = r{ai) 
and is defined by the local behavior of d near a; (see (21b)). 

After a conjugation of the jump matrix by diag{die ,dj^e '^"^), a fur- 
ther translation A \ — x/2 (to complete the square) and a rcscaling A 
wc end up with a problem with jumps on a rescalcd cross which is no more small. 
Extending the cross to infinity (by setting the jump equal to the identity on the 
extension) we have a new Riemann-Hilbert problem which can be approximated by 
the following: 

^+ = ^ -Vxs, where u*(A) 
Q ^ \ , for XeJi, 



(A.2) 



1 0\ . , - 

^-2u,i^2i\' I j , jor Xe J2, 



1 -nA^'^'^e-^'^' 
1 



1 

rtX-^^'^e^'^^ 1 



for X G J3, 

, for X& J4, 

limx^ooS{X) = I, 

where J, is the extension of J,, i = 1, 2, 3, 4. 
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The last Riemann-Hilbert problem can be explicilty solved via parabolic cylinder 
functions (see e.g. [DZl]). Indeed, let Da{X) denote the standard parabolic cylinder 
function. Then £>„(e-3''^/^A) and £)„(-e-3''^/*A) solve the ODE 



+ (o-^+«)^ = 0' 



rfA2 ^2 4 



where a = wi. Let 

(A.3) S = 

where 



Sii Si2 

S21 S22 



5n = e-3'''''/''Da(e-3"/4A), 
= (4lT.ll/4 ^"-^^[^^-^(^---^^A) - ^i^-„(e-/^A)], 

522 = e""'/^£'„(e-*^/^A), 
for ImX > 0, 

and 

5n = e'^'^'/^£)a(e''^/*A), 

52. = (,g|:;!,, e3--/2[Az^„(e^-/^A) + ^i^„(eW4,)], 

522 = e-^'^"'/*£'a(e3'"/^A), 
/or JmA < 0. 

Then it is possible to check that 

5+ = S^vs, where 5(A) = 

J forXGJu 



(A.6) ^ , 

1 -n 
1 



( ? ) , for Xe J2, 

\ l+|r,P ^/ 



/or A G J3, 



1 
rr 1 



for X e J4, 
limx^ooS{X) = I. 
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Now, setting 
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(A.7) S = ^X^^'^^e-'^^/^ 

and using the well known asymptotics for the parabolic cylinder function at infinity 
(A.8) 

D,{X) = A"e-^'/4(l +0(A-2)), \argX\ < ^, 
D„(A) = A«e-^^/4(l + 0(A-2)) - ^!ll^e"«A-«-ie^V4(i + o{X-')), \ < argX < ^, 

D,{X) = A«e-^^/^(l + 0(A-2)) - M^e-^-n— ie^^/^(l + 0{X-')), ^ < argX < 

it is immediate to check that 4* solves (A. 2). 

It can be easily verified that the back-rescaled local version near ai , 6; matches 
with the solution of the outer problem (as h ^ 0). Thus the issue of the small 
crosses is settled. 

APPENDIX 2. THE VARIATIONAL PROBLEM AND THE FINITE GENUS 
ANSATZ. 

The function g defined by (18) and the conditions stated before (18) is crucial for 
the asymptotic analysis of the Riemann-Hilbert problem (11). As stated in [KMM] 
and [KR] the existence of such a function follows from the existence and regularity of 
a solution to a variational problem. In this section we pose the variational problem 
and we state the results of [KR] on existence. We also show that a variation of the 
proofs of [KR] guarantees the validity of the finite genus ansatz for the barrier data 
problem. 

Let M — {z : Imz > 0}, be the complex upper- half plane and H = {z : Imz > 
0} U {oo} be the closure of H. Let also K = {z : Imz > 0} \ {z : Rez = 0, < 
Imz < A}, where A is a positive constant. In the closure of this space, K, we 
consider the points ix+ and ix-, where < a; < A as distinct. In other words, we 
cut a slit in the upper half-plane along the segment (0, iA) and distinguish between 
the two sides of the slit. The point infinity belongs to K, but not K. Define G{z; t]) 
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to be the Green's function for the upper half-plane 

and let dfi'^{ri) be the nonncgativc measure —p^drf on the segment [0, iA] oriented 
from to iA, where p° is the density of eigenvalues given by 

(9a) p°(A) = ^ 



7r(A2 + A2)i/2- 

The star denotes complex conjugation. Let the "external field" (j) be defined by 
(A.9) (l,{z) = - j G{z; ri)dijP{ri) - Re{i7v p°dr] + 2i{zx + zH)), 

where, without loss of generality x> Q. 

Let M be the set of all positive Borel measures on K, such that both the free 
energy 

E{ii) = j J G{x,y)dii{x)di^{y), ^gM 
and / (j)dfj. are finite. Also, let 

V^(z)= f G{z,x}dii{x), neM. 



be the Green's potential of the measure fj,. 
The weighted energy of the field (p is 



E^ilj.) = E{fi) + 2 j (f)dp., fi 



Now, given any curve F in K, the equilibrium measure supported in F is 
defined by 

E4,{X^) = min^^MiF)E4,{iJ.), 

where M{F) is the set of measures in M which are supported in F, provided such 
a measure exists. 

The finite gap ansatz is equivalent to the existence of a so-called S-curve joining 
the points 0+ and 0_ and lying entirely in IK. By S-curve we mean an oriented 
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curve F such that the equihbrium measure exists, its support consists of a finite 
union of analytic arcs and at any interior point of supp/j. the so called S-property 
is satisfied 

(A.IO) A.j,^v^'),A.j,^y^'y 
In [KR] we show that there is a C such that 

(A. 11) E^{\^) = max contour sfE^{\^) = maXcontoursFmin^^M{F)E(l,{lj) , 

and that the existence of an S-curve follows from the existence of a contour C 
maximizing the equilibrium measure. 

EXISTENCE THEOREM [KR]. For the external field given by (A.9), there 
exists a continuum F e F such that the equilibrium measure exists and 

(A.12) E^[F\{= E^{X^)) = maxFew'min^^M(F)E^{n). 

PROOF: See Theorem 4 in [KR]. Even though the density p does not satisfy all 
of conditions (1) of Theorem 4, the proof goes through unaltered. 
For our particular problem we also have 

REGULARITY THEOREM. Under the extra assumption that the continuum 
F does not touch the spike {z : Rez = 0, < Imz < A} except at a finite number 
of points, the continuum F is at worst a union of an S-curve and a finite union of 
real intervals. 

PROOF: The proof follows as in Theorem 8 of [KR]. But there are some changes 
here. The density of eigenvalues given by (9a) does not satisfy all conditions (1) 
set in [KR]. In particular it is not true that Im{p°{z)] > 0, for z G (0, iA] U IR+. 
Rather Im[p'^{z)] = on the real line. So a small amendment of the regularity 
proof is needed. 

The point of the assumption that Im[p°{z)] > 0, for z € (0,i^] U IR+ is to 
ensure that the continuum F does not touch the negative real line, except of course 
at 0_ and possibly infinity. In our case, we can argue as follows. 
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Because of the analyticity properties of the field, the real line can be divided 
into a finite number of intervals Jk such that in the interior of each Jfe either 

dimz ^ ^' 

or 

dImz ~ 

In the first case, one can see that the continuum F cannot touch the real line except 
of course at the endpoints of J^.. This is because for any configuration that involves 
a continuum including other points on the real line, we can find a configuration 
with no other points on the real line, by pushing measures up away from the real 
axis, which has greater (unweighted and weighted) energy. It is crucial here that if 
u e M then G{u, v) = 0, while if both u, v are off the real line G{u, v) > 0. 
In the second case, we get a finite union of real intervals. 

So F consists of a finite union of arcs: some of them are real intervals (at worst) 
and the rest do not touch the real line except at their endpoints. 

To pursue the proof of regularity one neeeds the following identity. 

THEOREM [KR]. Let F be the maximizing continuum of and be the equi- 
librium measure. Let /x be the extension of to the lower complex plane via 
H{z*) = -n{z). Then 

Rei[^ + V'{z))' = Re{V'{z))' - 2Re ( m^lMd^(^) 
J u — z J z - u 

+Re[j^ j 2{u + z)V'iu) dniu)]. 
Here V is the complexification of the real field (p. 

PROOF: By taking variations with respect to the equilibrium measure; see The- 
orem 5 of [KR]. Of course, variations cannot be taken near the real intervals, but 
since the field ^ = and the Green's function is also zero, these real intervals are 
not part of the support of the equilibrium measure. 

Prom (A. 13) it is easy to see that the support of the equilibrium measure of the 
maximizing continuum is characterized by 

(A.14) Re J {R^f^dz = 0, 



(A.13) 
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where 

R,iz) = iv'iz)r-2[ n^i^Mu) 

J suppfi ^ ^ 

+ 2{u + z)V'iu) diiiu)). 

Since R^iz) is a function analytic in K, the locus defined by (A. 14) is a union of 
arcs with endpoints at zeros of Rfj,. But it is readily seen that i?^ has finitely many 
zeros. 

(A. 10) also follows easily from (A. 13); see [KR]. Alternatively, see Chapter 8 of 
[KMM]. 

If F touches the spike [0, iA] at more than a finite number of points, regularity 
cannot be proved as above because variations cannot be taken. In [KR] we have 
included a rough idea on how to extend the above proof. Here is a more detailed 
argument forthcoming (see also [K]). 

One wishes to somehow allow the contour F go through the spike [0,iA]. One 
problem arising is that (the complexification of) the external field is not analytic 
across the segment [— iA, zA]. What is true, however, is that V is analytic in a 
Riemann surface consisting of infinitely many sheets, cut along the line segment 
[— i^jiA]. So, the appropriate underlying space for the (doubled up) variational 
problem should now be a non-compact Riemann surface, say L. Now, compactness 
is the crucial element in the proof of a maximizing continuum. But we can indeed 
compactify the Riemann surface L by mapping it to a subset of the complex plane 
and compactifying the complex plane. The other problem, of course, is whether 
the amended variational problem (with the modified field defined on the Riemann 
surface and with the possibility of F not enclosing all the original eigenvalues) is 
still appropriate for the scmiclassical NLS. The argument goes roughly as follows: 

(i) Proof of the existence of an S-curvc _F in L along the lines of [KR] . 

(ii) Deformation of the original discrete Riemann-Hilbert problem to the set F 
consisting of the projection of F to the complex plane. At first sight, it is clear that 
F may not encircle the spike [0,iA]. It is however possible to append S-loops (not 
necessarily with respect to the same branch of the external field) and end up with 
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a sum of S-loops, such that the amended F does encircle the spike [0,Zj4]. To see 
this, suppose there is an open interval, say {ia,iai), which lies in the exterior of 
F, while ia, iai e F. Let us assume, that F crosses [0, iA\ along bands at ia, ia\ 
(if not the situation is similar and simpler) ; call these bands S, S\ . Let (3~ , /3+ be 
points (considered in C) lying on S to the left and right of ia respectively, and at a 
small distance from ia. Similarly, let ,(3'^ be points lying on to the left and 
right of iai respectively, and at a small distance from iai. We will show that there 
exists a " gap" region including the preimages of (3~ , lying in the A''th sheet for 
—N large enough, and similarly there exists a "gap" region including the preimages 
of /?+ , (3i lying in the Mth sheet for M large enough, both being regions for which 
the gap inequalities hold a priori, irrespectively of the actual S-curve, depending 
only on the external field! 

Indeed, note that the quantity Re{(j)'^{z)) (which defines the variational inequali- 
ties) is a priori bounded above by —4>{z). For this, see (8.8) in Chapter 8 of [KMM]; 
there is actually a sign error; the right formula is 

Re{r{z)) = -<f>{z) + J G{z,rj)p'^{rj)drj. 

Next note that the difference of the values of the function Re{(l)'^(z)) in consecutive 
sheets is 6Re{<j)'^) = ±2-KRez, and hence the difference of the values at points on 
consecutive sheets whose image under the projection to the complex plane is ir] + e, 
where r] is real and e is a small (negative or positive) real, is 5{Re^'^) = ±27re. 
This means that on the left (respectively right) side of the imaginary semiaxis, the 
inequality Re{(f)'^{z)) < will be eventually (depending on the sheet) be valid at 
any given small distance to it. 

Applying the theory of [KR] wc join the preimages of /3~ and (3^ (under the 
projection pf L to C) lying in the A^th sheet and the preimages of /?+ and [3^ lying 
in the Mthe sheet. Finally we connect the preimage of /3+ lying in the TVth sheet 
to the preimage of (3~ lying in the Mth sheet, and so on. We thus end up with an 
S-loop whose projection is covering the "lacuna" {ia, iai). 

The original discrete Riemann-Hilbert problem can be trivially deformed to a 
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discrete Riemann-Hilbert on the resulting (projection of the) union of S-loops. All 
this is possible even in the case where F self-intersects. 

(iv) Deform the discrete Riemann-Hilbert problem to the continuous one with 
the right band/gap structure (on F; according to the projection of the equilibrium 
measure on F), which is then cxpUcitly solvable via theta functions. Both the 
discrete-to-continuous approximation and the opening of the lenses needed for this 
deformation are justified as in [KMM] (see also the article [LM] for the delicate 
study of the Riemann-Hilbert problem near the points where F crosses the spike). 
The g-function is defined by the same Thouless-type formula with respect to the 
equilibrium measure (cf. section 2(iii)). It satisfies the same conditions as in [KMM] 
(measure reality and variational inequality) on bands (where the branch of the field 
turns out to be irrelevant) and on gaps (where the inequalities are satisfied according 
to the branch of the external field). 

We arrive at the following conclusion. 

FINITE GAP ANSATZ. The finite genus ansatz is valid generically in x,t for 
the steplike data given by (2). 
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